Supplemental Table 1. Sample information 
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Table S2: Three-epoch demographic model results from dadi for individuals on either side of a 
putative cline. 

Inferences based on subsets of the western lowland gorilla samples, using the 8x data. Samples on the 
1 st side of cline: KB3782, KB3784, KB5792, KB5852, KB7973, X00109. Samples on the 2 nd side of 
cline: A930, A931, A933, A936, A937, A962, KB6039, X00108. 



Table S3: CLR Window 
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0.00825 



Table S4: Top enriched GO terms from TopGo. 

Results based on the Fisher's exact test and the elimination method by Alexa et al. (2006). All genes in 
regions with p < 10" 3 were included. The 16 terms with p < 0.01 are listed. 

Alexa A, Rahnenfuhrer J, Lengauer T. 2006. Improved scoring of functional groups from gene 
expression data by decorrelating GO graph structure. Bioinformatics. 22(13): 1600-1607. 
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Table S5: Comparison of published estimates of gorilla population split times. 

Values from each study are adjusted to match the generation time and mutation rates employed in this 
study. A mutation rate adjustment is not performed for Thalmann et al (201 1), which is based on 
microsatellites. We note that Thalmann et al (2007) report a range of split time values under several 
different models, with population split times ranging from 78 kya (without subsequent gene flow) to 
1.6 mya, with some models including Eastern to Western migration. 
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Figure SI: Genome wide heterozygosity values. 

Heterozygosity was calculated along the autosomes for each sample based on the 8x coverage masks. 
Samples are colored based on species: eastern lowland (red), Cross River (yellow), and western 
lowland (blue). 
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Figure S2: PCA of three gorilla species. Colors indicate species: Cross River (brown), western 
lowland (blue), and eastern lowland (red). Data was thinned to include 20% of the total number of 
intersecting SNPs. Percentages indicate the percent of variance explain by each principle component. 
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Figure S3: PC A of western lowland gorilla. 

PCA of 14 western lowland gorillas based on 8x sequencing data. Colors indicate geographical origin: 
Equatorial Guinea (orange), Cameroon (blue), Congo (red), captive born (black), unknown origin 
(grey). Percentages indicate the percent of variance explained by each principle component. 
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Figure S4: Admixture results for the three gorilla subspecies. 

A: Population membership inferred from ADMIXTURE with 2-4 populations. B: Cross Validation 
Error from ADMIXTURE analyses on three gorilla subspecies. CV error is estimated for 10 
independent runs, and the results from run with the CV error is shown in A. 
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Figure S5: Admixture assignment and CV Error for 14 western lowland gorillas. 

A: Population membership inferred from ADMIXTURE with 2-4 populations. B: Cross Validation 
Error from ADMIXTURE analyses on three gorilla subspecies. CV error is estimated for 10 
independent runs, and the results from the runs with the lowest error is shown in A. 
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Figure S6: Four-population phylogeny as inferred by TreeMix (Pickrell and Pritchard, 2012). 

The tree was inferred using reference-free genotypes obtained using BSNP at putatively unlinked 
neutrally evolving loci. 
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Figure S7: Admixture analysis for subset of samples used in G-PhoCS analysis. 

Results are based on reference-free genotypes obtained using BSNP at putatively unlinked neutrally 
evolving loci. 
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Figure S8: Total migration rate inferred by G-PhoCS, given by m AB * x AB . 

Along the x-axis values from different migration bands associated with different scenarios (Figure 1) are depicted. The first six migration rate 
estimates come from a bi-directional migration setting. The last setting is with two migration bands, one from western lowland to eastern 
lowland gorilla and another from Cross River to eastern lowland gorilla. A strong signal of migration of western to eastern lowland gorillas is 
observed. Colors correspond to different combinations of western gorilla samples included in the analysis. 




Figure S9: MCMC traces for three demographic parameters inferred using G-PhoCS. 

Each trace plot (obtained using Tracer) displays 9 combinations of samples including 2 eastern gorillas, 2 western gorillas and 1 Cross River 
gorilla with migration setting scenario 5 (Figure 1). A. Western gorilla population size B. Western-Cross River gorilla split time C. Eastern- 
Western-Cross River ancestral split time 
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Figure S10: Eastern lowland gorilla population size estimates by G-PhoCS. 

We include 2 eastern gorilla individuals (Victoria, 9732), different combinations of Western individuals, 1 Cross River gorilla (B646) and 
always include human reference as an outgroup. Values along the x-axis correspond to five different scenarios depicted in Figure 1 . Raw 
estimates are depicted on the left axis scale, and recalibrated values (effective population size) assuming a human-gorilla divergence time of 12 
Mya are depicted on the right axis. 
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Figure Sll: Western lowland gorilla population size estimates by G-PhoCS. 

Values along the x-axis correspond to five different scenarios depicted in Figure 1 . Raw estimates are depicted on the left axis scale, and 
recalibrated values (effective population size) assuming a human-gorilla divergence time of 12 Mya are depicted on the right axis. 
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Figure S12: Western lowland and Cross River gorilla split-time estimates by G-PhoCS. 

Values along the x-axis correspond to five different scenarios depicted in Figure 1 . Raw estimates are depicted on the left axis scale, and 
recalibrated values assuming a human-gorilla divergence time of 12 Mya are depicted on the right axis. 
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Figure S13: Eastern lowland and western-Cross River ancestral gorilla split-time estimates by G-PhoCS. 

Raw estimates are depicted on the left axis scale, and recalibrated values (effective population size) assuming a human-gorilla divergence time of 
12 Mya are depicted on the right axis. 
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Figure S14: Comparison of dadi and PSMC Models. 

The red line indicates population history obtained for a western lowland gorilla using PSMC based on 
mapping to the gorilla genome and scaled using generation time and mutation rates employed in this paper. 
The blue line corresponds to the best-fit three-epoch model inferyred for western lowland gorillas using 
dadi. 



